SURFACE SPLINE APPROXIMATION ON SO (3) 
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Abstr act. The purpose of this article is to introduce a new class of kernels on 
SO(3) for approximation and interpolation, and to estimate the approximation 
power of the associated spaces. The kernels we consider arise as linear com- 
binations of Green's functions of certain differential operators on the rotation 
group. They are conditionally positive definite and have a simple, closed-form 
expression, lending themselves to direct implementation via, e.g., interpola- 
tion or least-squares approximation. To gauge the approximation power of the 
underlying spaces, we introduce an approximation scheme providing precise 
L p error estimates for linear schemes, namely with L v approximation order 
conforming to the L p smoothness of the target function. 



1. Introduction 

Approximating a function / : X — > C on various structures X by linear com- 
binations of translates of a given single basis function is a widely used method. 
At the heart of this methodology is the creation of an approximant using a fixed 
kernel k : X x X — > C, by taking a scattered linear combination of copies of this 
kernel: ss(x) = X^es A^k(x, £). The success of these methods derive from their 
ability to generate approximants from data having arbitrary geometry. The case 
where the underlying structure are Euclidean spaces M. d or Euclidean spheres S d has 
been studied in great detail over the last decade, see [H [31] . Such kernel methods 
have found success in approximation theory (for treating high dimensional scat- 
tered data) and learning theory (in kernel based learning, support vector machines, 
neural networks, etc). 

However, in various applications we are confronted with the situation that the 
underlying set is a compact or locally compact group and one is asked to propose 
suitable approximation procedures on these structures (this is a setting with no 
natural dilation operator and no 'regular grids'). Such problems arise in biochem- 
istry, crystallography and robotics to name only a few. The monograph [4] provides 
a great collection of scattered data approximation problems where different matrix 
groups are involved. 

In view of applications, the rotation group SO(3) is, without doubt, one of the 
most important groups. There exist a wealth of scattered data approximation prob- 
lems on SO(3), cf. [27l EH [3j [32] ■ Such problems have attracted significant recent 
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attention from the mathematical community, which has investigated approximation 
via Wigner functions [13j O OH [HI [25] , and kernels PUCE]. 

The goal of this article is to introduce a family of kernels on the rotation group 
suitable for approximation and to investigate the approximation power of this family 
by estimating the performance of a (theoretical) approximation scheme. The kernels 
we present are notable in many respects: 

• they are conditionally positive definite, have a simple, closed form expres- 
sion, and are consequently well-suited for implementation (for, say, inter- 
polation or least squares approximation); 

• they invert iterated, perturbed Laplace - Beltrami operators, tying them 
directly to successful kernel based approximation in other settings (cf. sur- 
face spline approximation on R d jS], [7], restricted surface splines on spheres 
[15] , and, of course, polynomial splines in one dimension); 

• they have excellent scaling properties, much like the surface splines in the 
Euclidean case, suggesting that they can adjust themselves to data whose 
density varies spatially, as in [BJ. 

In addition to these three points, we study the approximation power of the spaces 
generated by the kernels. We develop an approximation scheme and accompanying 
error estimates derived from the kernel's role as a Green's function. To be sure, 
the scheme we develop is theoretical and does not immediately lend itself to direct 
implementation, but the accompanying error estimates are of a type that is notori- 
ously difficult for kernel based approximation when the kernel is not dilated. This 
is true even in the Euclidean case, as in radial basis function (RBF) approximation. 
The results we are after give precise error rates for functions at all appropriate levels 
of smoothness; we determine L p rates of convergence commensurate with the L p 
smoothness of the function being approximated. We note that the results we seek 
are well known for kernel approximation where dilations and regular grids are em- 
ployed [5j[2T], and even without grids, [20l [18]. However, in the absence of dilation, 
error estimates are often known only for target functions residing in a reproducing 
kernel Hilbert space. 

We have organized the article as follows. In the next section we give the neces- 
sary results of analysis on the rotation group. In Section 3 we introduce the kernels 
and give their expansions in terms of Wigner D-functions, the basic elements of 
Fourier analysis on the group. From this we are able to show that the kernels are 
conditionally positive definite, and, hence are capable of solving interpolation and 
other scattered data fitting problems. The expansion also allows us to identify the 
integral identities these kernels satisfy. Section 4 establishes the basic strategy to 
localize a kernel: by replacing a given kernel by a linear combination of its nearby 
copies. A precise estimate of the error in making this exchange is then given. In 
Section 5 we propose an approximation scheme and present error estimates that 
demonstrate the approximation power of the kernels. Finally, Section 6 involves a 
discussion of potential generalizations of this method to other groups and homoge- 
neous spaces. Recent work on bounding Lebesgue constants for kernel interpolation 
is discussed, which indicates that the theoretical error estimates may also hold for 
more practical schemes like interpolation and least squares approximation. 
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2. Background 

In this section, we collect some basic material and necessary notation on the 
rotation group to keep the paper self-contained. 

2.1. Basic Facts. Let 5*0(3) := {x G K 3x3 : x T x = e, deta; = 1} denote the non- 
Abelian compact group of proper rotations in the Euclidean space R 3 and let p 
denote the normalized Haar measure on 50(3), i.e., we have /i(SO(3)) = 1. There 
are various ways to parameterize this group. An important parametric form for the 
rotation group uses the so-called Euler angles ((pi,9,(p2) € [0,27r) x [0,7r] x [0,2-7r). 
In the literature there are miscellaneous conventions of Euler angles. In this paper 
we follow the conventions made in [ITJ Section 1.2]. There the Euler angles are 
defined such that we can write every x E 50(3) as 

X = x(tpi,6,Lp 2 ) = S z [ip\}3 x [6]8 z [ip2], (1) 

where 





are rotations about the z-axis and the x-axis, respectively. 

Any integral over 50(3) (or portions thereof) will always be assumed to be the 
Haar integral. Using Euler angles, the Haar integral of a function / on 50(3) reads 
as 

/ f(x)dp(x) = — / / f(<px,e,ip 2 ) sin0d<^d0d^ 2 . (2) 
Jso(3) 87r Jo Jo Jo 

Besides the well-known parameterization via Euler angles, the parameterization 
via the three-dimensional projective space is also of some importance to us. Let IC^ 
be the closed ball in R 3 of radius ir centered at the origin and identify antipodal 
points. This is the three-dimensional projective space. An element x € 50(3) is 
identified with a point in the projective space IC^ by x — > ljv where r, satisfying 
xv = r and ||r|| = 1, is the rotation axis and w, which can be chosen in [0, 7r], is the 
rotation angle of x. 

Then it is easy to see that 

dist(a;, y) := uj(y~ 1 x) (3) 

defines a translation invariant metric on 50(3). Of course, this is not the only met- 
ric one can place on 50(3), see [22] for a discussion of other metrics and sampling in 
50(3), but it is compatible with the Riemannian structure of the group and we will 
use this metric throughout this article. This leads to the usual definition of the ball 
centered at a having radius p: B(ct, p) = {C I dist(C, a ) < p)- A consequence of © 
is that we can estimate the volume of a ball p.(B(a, p)j = J so ^ ^-B(a.p) ( x )dp(x), 
obtaining constants b\ := 2/(37r 3 ) and &2 := 1/(6tt) so that for any < p < ir and 
any a € 50(3), 

bip 3 < p{B(a,p)) < b 2 p 3 . 
We utilize the notation ui in the context of the rotation angle and dist in the context 
of the metric. 

We can use this metric in order to quantify the distribution of the centers in the 

set 

E:={£, G5O(3):j = 0,...,M-l}. 
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To do so we define two parameters. The first one is the separation distance 

q = q{S) "o^Vi^^-y 

which measures 'clustering trends' in H C 50(3). On the other hand, the fill 
distance 

h = h(a) := max mm dist(£.,, y) 
yeSO(3) j=o,...,M-i 3 

describes the density of H in 50(3). 

2.2. Harmonic Analysis and Kernels. The basic building blocks in harmonic 
analysis on topological groups are the irreducible unitary representations. On 
50(3), it is possible to compute a complete set of irreducible unitary represen- 
tations, explicitly. More precisely, using the quasi-left regular representation of 
50(3) on L 2 {E> 2 ) one is able to compute for each £ e N Q a, (2£+l) x (2£ + l)-matrix 
D l such that these operators D , I € No, form a complete set of irreducible unitary 
representations of 50(3). The matrix elements Dt in the canonical basis are 
often called Wigner-D functions of degree £ and orders k and m. In terms of Euler 
angles these functions are given by 

Di m (x) = Di m (<p u 9^ 2 ) = e-^P£ m (cos0)e-™^, -t < k,m < I, (4) 

where the function Pj. m is given by 



pi 

k.m 



(t) = C(l - t )-^{i + t)-^^{(i -tf-\i + tf+ k ) 



with r - (- 1 ) Sm ' i {t-k)\(t+ m )\ 

Wltn O - 2 t {i _ k)l y (£+fc)!(*_m)!- 

We refer to [11] for all the details in these computations. By means of the 
Peter- Weyl Theorem the Wigner-D functions, the matrix elements of the irreducible 
representations, constitute an orthogonal basis of the L 2 (S'0(3)). Hence, every 
/ € L2(SO(3)) can be expanded in a SO(3) Fourier series 



f=J2 E V2TTTfi m D 



with 50(3) Fourier coefficients ff. m = \J2£ + 1 J so ^ f{x)D e km {x)d[i{x). On the 
other hand, to each eigenvalue, vi :— 1(1+ 1), of the Laplace -Beltrami operator A 
on SO(3), there corresponds an eigenspace of dimension (2£ + l) 2 with orthogonal 
basis (Dl m )k,m, i.e. we have for all £ G No and all k,m — —£, . . . , £ that 

ADl m =£(£ + l)Di jm . (5) 

We denote this eigenspace by He, and we refer to Tl n = 0™ =o Hi as the polynomials 
of degree no more than n. 

The kernels we study in this article possess a certain symmetry; they are func- 
tions of the distance between the two arguments, dist(cc,a) = u(a.~ 1 x) € [0,7r] 
only. By the bijectivity of sin (2) and cos (2) on [0,tt] we are entitled to express 

these kernels as k(x, a) = (f>(sm( ul ^ a 2 = ^(cos( "^° 2 where 0, -0 : [0, 1] — > 

!L There are benefits to both alternatives, as we shall see in the following sections. 

On 5*0(3) one easily checks that the functions that only depend on the rotational 
angle of the argument coincide exactly with the class functions - functions that are 
constant on conjugacy classes, i.e. f(x) = f(yxy~ 1 ) for all x,y e 5*0(3). In 



SURFACE SPLINE APPROXIMATION ON SO(3) 



5 



other words, for any class function / on 50(3) there is a uniquely determined 
/ : [0,?r] -> C such that f(x) = f(uj(x)) for all x e 50(3). The Haar integral for 
such a class function reads as 

f(x)d/i(x)= [ /(w(aO)dM(s) = - / f(t)sm 2 (l)dt. (6) 

SO (3) JSO(3) n JO \ Z / 

We make the convention that if / is a class function on 5*0(3) and if no confusion 
occurs, we subsequently drop the tilde. 

Since the characters q = TiD e , I £ No, of 5*0(3) form an orthonormal basis for 
the space of class functions in 2^(50(3)), we can decompose any continuous kernel 
possessing the symmetry mentioned above as k(x, ex) — YlTLo k(£)ce(oc x). Using 

the formula Ce(at~ x x) = ^2£(cos( "^" 2 x ^ )), such kernels can be expressed as an 
expansion of even powered Chebyshev polynomials of the second kind, 



oo 



k(x, a) = k(£)c e (a- L x) = }_J k(£)U 2l (cos X) ) ) . (7) 

t=o i=o \ \ z J / 

Alternatively, we can write: 

oo oo 

k(x,a) = ^h^ceia^x) = ^rk(£)Tr(D e (a.- 1 x)) 

1=0 i=0 
oo 

e=o 

oo e t 

- E E E (8) 

e=o l=-i v=-t 

This is sometimes called the addition formula for Wigner D-functions. 

The expression of the kernel as a series of Chebyshev polynomials (J7]) as well as 
the addition formula (j8J) are of prime importance in the study of kernels on 50(3). 
The first is used to understand the kernel as the fundamental solution of a simple 
differential operator, which is the goal of the next section, while the second is key 
to localizing the kernel, which is tackled in Section 4. 

2.3. Smoothness Spaces. Because we are interested in approximating in the 
L p (SO(3)) metric, we must introduce L p smoothness spaces on 50(3). This is 
done in two stages. In the first stage, we introduce the (simpler) Sobolev and 
Holder spaces, whose definitions are accessible, if not familiar to most readers. The 
second stage is postponed until Section 5, where we introduce Besov spaces. 

Sobolev spaces on manifolds have a variety of equivalent characterizations. Per- 
haps the most straightforward of these is to import the definition from Euclidean 
space via a partition of unity and a corresponding set of diffeomorphisms. This is 
the approach taken in [26l 1.11]. The Sobolev space on a region O of E 3 is defined 
to be the space of functions / g L v {0) such that the quantity 

Wf\\w }{ o) : =E / l^/(aOl p d* (9) 



is finite. To transport this definition, let (V'j)jLi be a partition of unity of 50(3), 
and let (flj, hj)^L± be a corresponding collection of charts (each fij is an open set 
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in SO (3) containing the closure of supp(V>j) and each hj : ilj — > Oj C R 3 is a 
diffeomorphism) . 

Definition 1. For 1 < p < oo the Sobolev space W k (SO(3)) consists of functions 
f G L p such that 

ll/C* E ii to ° (hj 1 )) (/ o (hf)) \\ p Wpk(0j) 

3=1 

is finite. 

In this way we transport the Sobolev space definition from Euclidean space. For 
smoothness measured in L^, although the definition makes sense, we end up using 
the C k spaces, whose norm is defined according to the same procedure. 

An important observation, of use in future sections, is the mapping property of 
the Laplace - Beltrami operator. For j < k/2, 

II A'/IIp < ll/llwj. 

This follows from the well known expression of A in local coordinates [26] 7.2.5]: 
namely (A/) o h = L(f o h) where 

is a second order differential operator having smooth coefficients (here g is an 
invertible matrix having smooth entries and g-* k is the j, k entry of its inverse). 

Another useful observation concerns the density of C°° functions in Wp(SO(3)), 
(and in C (SO(3)) when p — oo), which follows from the density of such functions 
in K 3 . Each real valued function (V'j/) ° {hj ) can be approximated with error e in 
W k by a function gj : Oj — »• M (via mollification, for instance). The desired smooth 
function is 

N N 

< : ■■ E & - g j : = E fa x (ft ° h j) ■ 

3=1 3=1 

Here we utilize a second family of smooth functions: each ipj : 50(3) — > R is a 
C°° cut-off function satisfying supp(V'j) C Slj and i>j(ot) = 1 for a e supp(-0j) 
(the function ipj x Gj is extended by outside of Slj). It follows that G — f = 

T,jLi$j( G j -4>jf) and 

\\G - f\\w£{SO(3)) 
N N 

- EE I P X ^ X {ft ~ ^3 X /) ° (V 1 )} ° ^3 

*=1 3 = 1 
N 

< E ^ 

«,i=i 

where the last line follows from the fact that each hj o hj 1 is a diffeomorphism 
from hg(flj n n^) to hj(Slj (1 fi^). Therefore, ||G — /II w p fc (SO(3) < C £ > where C 
depends only on the two families of smooth functions (V'i)jLi an d (ipj)j-i, and on 
the diffeomorphisms fe,-. 




V>* x ipj o h 



c k <n,') 



\9j~to x/m^- 1 )! 



WHOA 
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3. Surface Splines on SO (3) 

In this article we investigate approximation properties of the surface splines 
(perhaps more accurately the 'rotational surface splines'), which are the functions 
derived from <f> s (t) :— t 2s or, alternatively tp s (t) = (1 — t 2 ) s with the index, s, 
a positive, pure half integer (which gives the order of zero of the kernel or its 
smoothness) related to the 'order' of the kernel by s = s(m) = 2 ™~ 3 . We have, 

M*,a) := (™( t ^)) ! "'" = *.(-( i ^)) (10) 

To ensure continuity of the kernel, we will assume throughout the article that 
to > 2. 

An analysis of their Wigner D-function expansions will reveal two remarkable 
qualities. First, the kernels are conditionally positive definite; this property and its 
consequences are discussed in 13.21 Second, they satisfy integral identities, valid for 
/ <= W* m (SO(3)), of the form 

/ = / C m f{a)k m (-.a)dfi(a) (11) 

JSO{3) 

where C m is an elliptic differential operator of order 2to on 50(3). To be sure, there 
are many kernels satisfying a reproduction formula like and many of them are 
conditionally positive definite. However, in general, it is unclear that they will have 
a 'nice,' closed form expression similar to (|10[) . Similarly, there is no guarantee that 
they will have the property of localizability (developed in the next section), which 
is key to our understanding of the approximation power of the kernel (as discussed 
in Section [5}. 

We organize the section as follows. In 13.11 we expand (fT0|) in terms of Chebyshev 
polynomials. This expansion is used to show that the kernels are conditionally pos- 
itive definite, which leads to a discussion of some practical approximation schemes. 
This is treated in !3.2l In l3.3l we establish the reproduction formula (fTTj) and demon- 
strate that the underlying differential operator C m is of the form p(A), a degree to 
polynomial in the Laplace-Beltrami operator. 



3.1. Chebyshev coefficients of the surface splines. Our investigation of k m 
begins with studying its decomposition in characters (equivalently, the decomposi- 
tion of ip s in even Chebyshev polynomials) . 

Lemma 2. For m > 2, the £ th coefficient of k m in its expansion in terms of 
characters, cf. {7]), is 

~^)=It^M n p+i+r 1 - (12) 
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Proof. We consider k m (x,ot) — k m (a. 1 x) — J2eLo^"i(£)ce(a 1 x) and we want 
to determine the coefficients k m (£). Utilizing the fact that the characters are or- 
thonormal we wish to compute 



k m (£) = / k m (x)c e {x)dfi(x). 

JSO{3) 

By a change of variable, we get 

Ufl-§/>(*))~^«(2)M>. 

Rewriting U^t (cos (^)) in the familiar form 

'w\\ _ sin ({21+ 1)|) 



^( C0S UJ) = sin (|) 



and applying another change of variable 6 — lo/2, we obtain: 



4 



tt/2 



k m (£) = - / (sin(9) 2m ~ 2 sin((2£+ l)0)d0. 

Such integrals are easily attacked using simple trigonometric identities. They can 
also be reduced to the observation that repeatedly multiplying a periodic function 
by sin# is equivalent to applying a difference operator to its Fourier coefficients. 
Doing so, we obtain 

( S in0) 2M S in(L0)= (4) M f>l)' ( 2 ^)sin([L-(2Af- 2j )]0). 

Applying this identity with M = m — 1 and L = 21 + 1, and integrating each term 
sin([2i? + 3 — (2m — 2j)]6) over the interval [0, we see that 

~ _ / y-Uy (-!)> /2m- 2 

mW ~ I 4j tt ^ 2£ + 3-(2m-2,)l J 

1 \ 2 2 ^ 2 (-l)J /2m - 2 

Hence, k m (£) is a difference operator applied to the rational function x i-> a; -1 . We 
simplify this type of expression in the elementary formula: 

yv_iy ( M ) _J_ = ^ da) 

^ Vi/i+i (i)(L + l)...(i + M-l)(L + Jlf)" 1 ' 

Thus, we can simplify k m (£), by using (IT3|) with L = £ — m + | and M = 2m — 2: 

* <£)-*( T" 1 (2m " 2)! 

mW 7T v 4/ 771 + 3/2)... (* + m- 1/2)' 

which proves the lemma. □ 
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3.2. Treating data with surface splines. At this point we observe that the 
kernels k m are each conditionally positive definite of order £q > m — 2, and, hence, 
well suited for solving interpolation problems. Being conditionally positive definite 
of order £q means that for any H C SO(3), X^e^ X^e^ a €^C^€C > f° r an y 
a ^ in the subspace defined by Ba = 0, where 

A:= [M€,C)]e,<eB and B ■■= l D lM]e=o,...,e ,-i<L,v<e;Ses- (14) 

Indeed, conditional positive definiteness follows directly from the fact that k m (£) 
does not change sign for £ > m — 2 (see [§J Proposition 4.3]). 

We now present three classical, practical approximation schemes involving the 
kernel k m . Namely, we approximate by functions of the form 

la I 

Ces f=o i,u=—t 

where the kernel coefficients satisfy 

^a^„({) = 0, 0<e<£ ;-£< i,v<L (16) 

Interpolation We define the interpolant to a function / at the centers S as the 
function of the form (|15[) where the interpolation conditions 

a/,B«) = /(*), ^S, 

and auxiliary conditions ()16j) are satisfied. 

The interpolant is determined by solving the system 

A B T \( a\ ( f 
B J I /3 J 1 



(17) 



where f = (/(£))^es- It is well-known that system (|17j) is uniquely solvable for any 
set of centers B that is unisolvent for the polynomial space II^ D , see [31] Chapter 
8.5] (we provide sufficient conditions for such sets, in terms of the fill distance, in the 
next section). Moreover, for our kernels k m , setting up this system is particularly 
easy (see Note Q] below) . We refer the interested reader to [311 E] f° r background 
on interpolation by conditionally positive definite kernels. 

Tikhonov regularization For treating noisy data, a common approach is to find 
the kernel approximant minimizing a certain quadratic form. For noisy data (y^)^eB 
sampled at S C 5*0(3), we wish to solve the problem 

mm I K£) - yd 2 + M s \% 




W = T,t>t T,-e< h u<t \s:,u\ 2 \k m (i)\ 1 and a smooth- 
ing parameter A > is chosen based on the level of noise inherent in the data. The 
minimizer is of the form (|15|) . Because k m is conditionally positive definite, the 
minimizer is determined by solving the system of equations 



A + Aid B T \ f a \ _ ( y 



B J \ (3 J \ 

See [30l 2.4] for a detailed discussion of this problem. 



(18) 
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Least squares approximation For a unisolvent set S, the space 



S~ 



a^km(-,$,) + P | P e He , {tt$)$es satisfying {H 



has dimension #3, and one can construct the least squares operator P=, which is 
the orthogonal projector P& : L 2 (SO(3j) -+ 5= with range 5s- 

One may calculate the least squares projector by using any basis (vg)^ for 5s (say 
the basis of Lagrange functions \£ where xs, (C) — <^£,c)- ^ n ^ ms case J the projector is 
simply P=f = Y^t^B a £V£, with coefficients a = (a^ggs given by a = G -1 /, where 
G%,C = ( v £i v c) is the Gram matrix for the basis and and / = Since 
the elements of Ss are continuous and 50(3) is compact, the projector extends 
as an operator to each space L p (50(3)), 1 < p < 00. The L p stability of such 
operators on compact Riemannian manifolds has recently been investigated in [16] ; 
an open problem, |16l Conjecture 3], is whether the stability and approximation 
results of that article hold for the kernels presented here. 

Note 1. Suppose we are given centers H = {£ • e 50(3), j — 0, ...,M — 1} in 
terms of Euler angles (as described in Section l2?Tj) . i.e., each £ ■ is determined by 
(t^ij, 0j, f2,j)- It is easy to show that we can write the distance of two centers ^ i 
and £ ■ in terms of Euler angles; in fact we have 



cos 



<Pl,i - <fl,j 



cos 



V 3 2,j - <f2,j 



cos 



, fl,i - <fl,j \ . ( f2,i - ¥2, 3 , 

sm sin — ■ — cos 



Thus, we can directly write down the closed form expression for the entry Aij of 
the matrix A as 



Furthermore, the nonequispaced Fourier matrix B can be computed directly from 
(HJ) . Such matrices have been studied in great detail and can be handled efficiently, 
see [13] . 

On the other hand, if the centers are given via their rotation axis and rotation 
angle, we can immediately compute the corresponding Euler angles, using the direct 
relation between the Euler angle and the axis- angle parameterization of 50(3), see 
[28l Chapter 1.4.4]). By doing so, we can again easily set up the matrices A and 
B as shown above. 

3.3. Integral representation. Returning to our investigation of the coefficients 
k m (£), we observe that k m (€) can be simplified by arranging factors symmetrically 
around £ + 1/2. Doing this, we have: 



k m {£) 



2 (2m- 2)! 

7T (-4)™- 1 



n 



1 



L (£+1/2) -j 

m — 1 



1 



1/2 



n 



1 



L [l + 1/2) +3 



2 (2m -2)! 1 

~ (_4)m-i (^+1/2) JLJL (£ + 1/2) 2 -f 



n 



(19) 
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The last line comes from rearranging terms, and noting that [J — k] [ J + k] = 
J 2 — k 2 . Since (£ + 1/2) 2 = £(£ + 1) + 1/4, we can rewrite the expression in terms 
of 1(1 +1) and obtain (after reindexing) 

^ - ^f^^ +i )n ^ +1) -V-i/4] - (20) 

The coefficients k m (£) involve the symbol for A, cr(A) = £(£+1). In fact, these are 
the reciprocals of a polynomial in the symbol, p(cr(A)) , which is nonvanishing, 
because the numbers [j 2 — 1 /4] never coincide with eigenvalues of A. So the operator 
induced by k m inverts the elliptic operator p(A). This is the point of the following 
lemma: 

Lemma 3. The kernel (x, a) <-> k m (x, a) satisfies the formula 
f = k m (-,a) C m f{a)dfj,(a) 

JSO(3) 



using the operator of order 2m, C m := ^2^-2)! Iljlo 1 ^ — r i) w ^ r j := J 2 ~ 1/4; 
and which holds for all f e W 2m (SO(3)) (or C 2m (SO(3)) when p = oo). 

Proof. For a sufficiently smooth / (say C°° initially), we have the absolutely con- 
vergent series / = YltYli v^iv^l.v Using the addition formula (|SJ), property (O 
and the orthogonality of the Wigner D-functions we get 

k m (x, a.)C m f(a) d/z(a) 

SO(3) 



/ y^,k m (£)ce(a 1 x) x 

JSO(3) t 



7T(-4) 



m-1 m - 1 



II (A- rj )x;E^^( a ) d M«) 



(2m- 2)! 

v ' j=0 (.',1/' 

7sO(3) f 



m-1 '"-i 



EE4'4^W II (^(^ + l)-r,)< 1 ,(«)dM«) 

£' (.',1/' ^ j=0 

~ ttC— 4.~\™—l m_1 1 

= EE^ffl ( 2m L 2V IT [^+i)-^]^-r T <^N = /W- 

The reproducing formula follows for all / by a limiting argument, using the density 
of C°° in Wp ln for 1 < p < 00 (or C 2m when p = 00) and the boundedness of 

Cm. ' □ 



4. Localizing the Kernel 



For a set of centers S C SO(3) we now wish to investigate a 'coefficient kernel' 
a : S x 50(3) — > R that will allow us to easily replace k m (x, a) with k(x,a) :— 
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X^£ ea et)k m (x, £) in the representation To do so, we must estimate the 

cost of replacing the kernel, ek m , given by the 'error kernel': 

e km (x,a) := \k m (x,a) - a(£, a)fc m (jc,£)|. 

€63 

This is a bounded, rapidly decaying function, inducing an operator on L p with 
small norm. 

The outline of this section is as follows. In 14.11 we develop conditions on the 
coefficient kernel a(£, a) necessary for proving our results, and it is shown that 
such conditions are satisfied for sufficiently dense centers. A practical procedure 
for producing the coefficients is given in 14.21 In 14.31 we estimate the decay of the 
error kernel ek m (x,a). 

4.1. Coefficient kernel conditions. The two key quantities we need to resolve 
are the polynomial precision (the degree of Wigner D-functions reproduced by the 
coefficients a(-, a.)) and the rate of decay of the error kernel. As in the Euclidean 
and spherical setting, these are related: the greater the polynomial precision, the 
more local the replacement error. 

Definition 4 (Coefficient Kernel Conditions). For a set of centers H C 50(3) 
there is a number p = p(H). reflecting the density of a in 50(3), and a measurable 
kernel a : H x 5*0(3) — > K satisfying the following three conditions: 

CKC 1 (Support): a(£, a) = /or dist(a, £) > p. 

CKC 2 (Precision): For any polynomial of degree at most L, i.e., for any 
p £ 11^, the following holds: 

y } j a(£, ot)p(£) = p(a) for every a £ 50(3). 

CKC 3 (Stability): There is a constant K > 1 such that for every a G 
50(3), we have 

X>(€>eOI <K. 

€63 

As in the Euclidean case, such a local polynomial reproduction property holds 
for sufficiently dense centers. In order to show this we need the following fact about 
polynomials on 50(3). 

Proposition 5. Let p € 11^. Then p restricted to a geodesic segment on 50(3) 
gives a trigonometric polynomial of degree at most L. 

L i 

Proof. Let p = J] a l V D[ „ £ Hl and an arbitrary geodesic segment G — 

{xqS v (lu) : oj e [0,t]} of length t € [0, 2tt) on 50(3) be given. Here s^(w) denotes 
a rotation about the rotation axis rj € § 2 with rotation angle u£ [0, 2ir). Then we 
can find y £ 50(3) such that s v (uj) = y~ 1 s z (uj)y for all uj e [0,t]. The rotation 
s z (uj) can, for example, be represented by the Euler angles {w, 0, 0). Thus, we get, 
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by using the representation flU, for x(ui) — xqs v (uj) £ G 

l e 

p(x(ui)) = p(x y~ 1 8 z (u))y) = ^2 E ^-^(^oST^My) 



TriAjD^xoy-'s^y))) 



L 

- E 

L 

= ^2Tr(D e (y)AjD i (x y- 1 )D e (s z (u;)) 

Tj 

= EE «^H) = EE^- ,lu 

i=o l,v=-i e=o L=-e 

We note that the operator Bi and its entries bf v are independent of w. The final 
equality is a consequence of <(4j) . □ 

The Coefficient Kernel Conditions are easily satisfied. The following lemma 
shows that for sufficiently dense centers, one can find a suitable kernel. 

Lemma 6. Given L, an integer reflecting the desired precision, and given centers 
H having fill distance h < ho (where ho := -^j^, a constant depending only on L), 
there exists a coefficient kernel a : S x 5*0(3) — » R satisfying the GKCs with radius 
p = CL 2 h (with C > a global constant) and K = 2. 

Before providing the proof, we remark that the restriction p < tt forces h < jjj-^ , 
which is the only restriction placed on h. 

Proof. Constructing a coefficient kernel follows a 'norming set argument', which 
relies on being able to capture the norm of a polynomial over a ball by its restriction 
to the centers inside the ball: 

ll P IUoc(B(p Q ,ri) ^ 2 H P | S n B (p ,p)ll^c(HnB(p ,p))- (21) 

Utilizing this to obtain a local polynomial reproduction is a fairly standard pro- 
cedure, and the basics of our argument can be found in [T51 Lemma 4.2] or [32 
Chapter 3]. So all we need to show is ([2"Tjl . 

The inequality (l2~Tj) relies on locating a center £ £ S near the point where the 
uniform norm is attained, p £ B(p ,p), so that a geodesic segment of length ui, 
starting at p and passing through £ lies entirely within B(p ,p). Specifically, we 
require 

distM 1 

oj ~ 16L 2 V ' 

The geodesic segment is determined by the curve 7 : [—u/2,w/2] H> S0(3), param- 
eterized by arclength. The main estimate we employ is Videnskh's inequality |29j 
which reads as follows. Let T be a trigonometric polynomial of degree at most L, 
then 

\T'(t)\ < 2L 2 cot(cj/4)||T|| ioo( _ w/2ia)/2) for w/2 < tt, |i| < w/2. 

The trigonometric polynomial t = F o 7 is the restriction of P, the polynomial 
in question, to the geodesic segment. We then observe that t(—u>/2) — t(—uj/2 + 
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dist(p,£)) < J_^ / / 2 2+diSt(P:€) \T'(t)\dt. Applying Videnskii's inequality 

Pip) ~ P<£) < di S t(p,0||r|| Loo( _ /2 , w/2) < -P(p). 

This gives us the desired inequality: P(p) < 2P(£). 

The fact that ([22jl holds follows by employing a cone condition. On 50(3), a 
cone c(q, 9, is, r) is the union of geodesic segments of length r emanating from q 
and having the initial tangent vector within 9 of the tangent vector v. We call 
9 the aperture and v the direction of the cone. It is an elementary fact that 
every ball in 5*0(3) satisfies a uniform cone condition. I.e., there is 9 > so that 
for any p Q 6 5*0(3) and p > and for every point in the ball B(p ,p), the cone 
c(p, 9, is, p/2) C B(p , p), where v is the tangent vector at p of the geodesic segment 
connecting p and p . 

An appeal to geometry shows that there is a global constant C depending only 
on 9 and the geometry of 50(3) so that a smaller ball b C 50(3) having small radius 
h can be placed entirely within the cone c(p, 9, is, p/2) with its center a distance of 
C'h from p. Thus there is a center £ that is a distance Ch < (C + l)h from p, 
for which the geodesic segment connecting p and £ of length p/2 is entirely within 
B(p , p). Therefore, p = 32L 2 Ch suffices. □ 

4.2. Constructing coefficients in practice. Although the previous lemma guar- 
antees the existence of a coefficient kernel sufficient for our purposes, it does not 
provide a viable method for generating such kernels in practice. To do so, we make 
the extra assumption that the centers are quasiuniform: that h and q, as defined 
in Section 2, are related by h/q < m for a fixed 'mesh ratio' m < oo. Armed with 
this assumption, we can estimate the number of centers in a given ball of radius 
p = CL 2 h (as was used in the previous lemma): 

# 3a = # ( p(a, P ) n s) < Mg^l! < g) 2 CW. 

We may determine a modified coefficient kernel, which we call a that satisfies the 
coefficient kernel conditions by taking the t-i minimizing solution at each a of the 
system Ba = D a , where {^> a )(j,k,e) = D%{ a ) is the vector of Wigner D-functions 
evaluated at a, the matrix B is the 'non-equispaced Fourier matrix' introduced in 
(TTJ}. I.e., the matrix whose rows are the D-functions evaluated at the various £ in 
S„. I.e., MU = D%{£). 

It follows that the £ 2 minimizing solution is given by B* (BB*)~ 1 D a , so we 
obtain a coefficient kernel a(£, a) satisfying conditions CKC 1 and CKC 2 of Def- 
inition|4]via (a(f,a)) t „ = B* (BB*)~ 1 D a and zero extension. Moreover, it is 
clear that on the subset 

K e 50(3) | S„/ = 3 a } 

(i.e., the set of a' for which the sets S a / coincide) each a(£, •) is a Wigner D- 
function. Hence, a is measurable. 

To show that a satisfies CKC 3 we can compare it to the kernel a guaranteed by 
Lemma [S] to observe that 

||S(-,a)|U l(a) < (#H a ) 1/2 ||S(,a)|| £2(3) < V (tt/2)^L^ ||a(-, a)|| /a(s) 
< ^VC 3 L e m 3 ||a(-, a)||^ ( s) < nVC 3 L e m 3 . 
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So CKC3 is satisfied with K = WC 3 I 6 m 3 . 

4.3. Estimating the cost of replacing the kernels. A consequence of the Co- 
efficient Kernel Conditions is that for any x e 5*0 (3) and a smooth class function 
k, there is a convenient mechanism to estimate the error e K (x,a) in terms of an 
interval around uj{ot.~ l x). For fixed x G SO(3) we choose to express the kernel k in 

terms of t x (£) = cos( "^ 2 x ^ ), by way of the formula k(x, £) = ^(tJ(C)) assuming 
$ is smooth on the interval I x := [minQxjinaxQ,,,], where the set Q x depends on 
the rotation angles used, Q x :— {t x (ol)} U {t x ($,) : £ e (H n B(a, p)}. The cost of 
replacing the kernel can be estimated in terms of the length of the interval X x and 
the size of derivatives of ■& purely on I x . This is the gist of the following lemma: 

Lemma 7. Given a kernel k(x,C) = ^(t^(C)), with i9 € C L+1 (1 X ), and a coeffi- 
cient kernel satisfying the CKCs with precision L, the replacement error e K (x, a) = 
\K(u(a~ 1 x)) — X^e- a (£> a)K(u}(£~ 1 x))\. We have the estimate: 

e K (x, a ) < \U L+1 maxl^WI. (23) 

Proof. Let both x and a be fixed, and choose the Taylor polynomial of degree L, 
QL.t^ia), °f & expanded about t x (a) — t x (ol). Now qL.t m (a) ma Y be rewritten as a 
linear combination of even degree Chebyshev polynomials, 

L 

qL, tas ( a ){t) = J2c e U 2 e(Vt). 

1=0 

This follows because Uit is of exact degree 2£ and is a linear combination of even 
degree monomials only, so Un{\fi) is a polynomial of exact degree £, and, thus, the 

set (bt2e(Vt)) £ _ is a basis for set of polynomials up to degree L. Note, furthermore, 
that qL,t a (a)(tx{ci)) - 52a-(£,<x)qL,t*(a)(tx(£)) = by the addition formula, since 

L 

£=0 1<L k,m 

and each D l k m (C) is annihilated by fi = S a — a (£i a )$e- 

Consequently, e K (x,a) < £ € l a (£, a)| \ti(t x (£))-q Li t a ( a )(tx(€))\, and the lemma 
follows from the remainder in Taylor's theorem: 

W(t x (Z)) -gL, t . ( a)(t.(0)l < jLTYy\t*(t) -t.(a)| i+ W|0< i+1 >(t)|. 

□ 

We now measure ek m , the cost of replacing k m (x, a.) by a linear combination of 
its copies, a(£, a)fc m (a;,£). 

Lemma 8. Assume a is a coefficient kernel satisfying the CKCs with radius p, 
precision L > 2m and stability constant K = sup Q ||a(-, a)||^ 1 (s)- Then the re- 
placement error of the kernel k m satisfies 

e km (x, a) < C(L, m) Kp 2m ~ 3 (l + dlst( J' a) ) ^ * \ 
w/iere C{L,m) is a constant that only depends on L and m. 
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Proof. We break this into two parts, the first concerns x near to a (say dist(x, a) < 
2p) where we can make use of the high order zero of the kernel, while the second 
treats x sufficiently far from a so that we can utilize the smoothness of the kernel 
and Lemma [Jj 

Case I: Assume dist(cc,a:) < 2p. This implies that dist(a;,£) < 3p for every 

and using 



dist(:c,£) 



| G supp(a(-, a)). Writing the kernel as k m (x,£) = [sin 2 
the stability of the coefficient kernel, we see that 



\k m (x,a) -^2a(£,a)k m (x,£)\ < (l + iT)(sin-p 



2m -3 



2m-3 



< 



(1 + JQ(- 



2m -3 



„2m-3 



Because 1+K < 2K < (3/2) 2 K, it follows that e km (x, a) is less than {3/2) 2m ~ 1 Kp 2m - z . 
The condition that dist(cc, a.) < 2p implies that 3 > l + dist(a;, ot)/p, which handles 
the first case. 

Case II: We use the fact, from the definition of the surface spline (fT0|) . that 
k(x, £) = (!- cos 2 dist(c ' a;) ) s with 2s = 2m - 3. Since < dist(C, x) < it, we can 



rewrite the kernel as k(x,£) = $ s (tx(G))> with t x (C) ■= cos 2 ( 



2 / dist^aj) \ _ 



.(0)' 



and $ s := (f — -) s . We note that fl s is smooth for t less than 1. Lemma [7] tells us 
that the error is controlled by |Ia;| i+1 . The length of X x can be written as 



!*«(€) -*« (01 



cos 2 (^1^1 

. ( distff, x) 
sm - 



— sin 



2 ( dist(C, x) 

■s 

V 2 

/ dist(£, x) 

V 2 



. /dist(|^)\ . ^dist(C,g 
m + sm 



Thus, for x in a sufficiently wide annulus about a, say 2 n p < dist(a:, a) < 2 n+1 p, 
we have 2 n ~ 1 p < dist(cc,£) < 2 n+2 p for all £ G supp(a(-, a)), and we can estimate 

%,\ by 

|dist(C,OI |dist(as,C) + dist(x,C)| 



Applying Lemma [7j gives 



< P x (2»+ 2 p) 



e km (x,a) < 



K 



(L + l)V 



(L+l 



Differentiating gives \\$ 



lUoo(Xx) S 



) IU 00 (i„)P 2(i+1) 2 ( " +2)(L+1) . 



n 7 =o( s_ i) max teio,( 1 -*) s \ and we 



cos 2 dlht (£' a! ) for the rotation 



note that the maximum is attained at t = ^(C) 

C € {a} Usupp(a(-, a)) nearest to x. Rewriting f — t x (£) as sin 2 dlst( £> x > | anc j U ging 
the inequality sinfl > §0 gives f - t x (C) > (f dist ^> ) 2 > i^^) 2 - Thus, 

s-L-l 



/ 9 2(n-lJ n 2\ 

|^ +1) IU.(x,) < |a(— 



| S (S - 1) — (S - L)| ^2(L+l- fl ) p 2( S -L-l) 2 2(„-l) (l 



-L-l) 
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Together with the estimate of the length of \I X | , we obtain 

e km (x, a ) < K i8 ^ + ^ S) \s(s - 1) - (a - L)\p 2 ("+^*-^). 
Because 2"+ 2 > 1 + dist( p x -° ) , the result follows. □ 



5. Approximation Scheme and Main Results 

We are now ready to introduce the approximation scheme central to our error 
estimates. Given a coefficient kernel a satisfying the CKCs with precision L > 2m, 
we define the operator T H on W 2m (SO{3)) by 

T s f(x) = J2 A tkrn(x,£) (24) 

with A% :— f SO ( 3 -) 'Cr?x/(«)a(C a)dfj,(ot) (recalling that the operator L m was defined 
in Lemma [3]). 

The approximant 

T s f:= f y)a(€,a)fc TO (-,€)W(«)dM(a) 

JSO(3) 

is the function obtained by replacing the kernel in the integral representation by 
linear combinations of its nearby shifts. 

An apparent drawback of this approach is the requirement that / should be in 
the domain of C m : this is a restrictive smoothness condition on the target function. 
To remedy this, we must modify the approach slightly, by splitting the (less smooth) 
target function into two parts: a component of sufficient smoothness (in L p ) and 
a rough but negligible component. The approximant to a general function will be 
obtained by applying the operator T= to its smooth part only. Precisely how this 
split is accomplished will be set aside for now, and revisited later in this section. 

5.1. Functions of Full Smoothness. We obtain error estimates by bounding the 
norm of the operator induced by the error kernel, i.e., 

E km : L p (SO(3)) L p (SO(3)), E km g{x) = [ e km (x, a)g{a)dfi(a). 

JSO(3) 

The bound on such an operator controls the error of the global approximant, as we 
will see in the following theorem. 

Theorem 9. Let m > 2, and let a : S x SO (3) —> K be a coefficient kernel 
satisfying CKCs with radius p, precision L > 2m, and stability constant K . For 
f £ W p m (SO(3)), 1 < p < oo, the approximant converges, in L p , to f with error 
bounded by 

\\f-Tef\\ p <Cp 2m \\f\\ w}m{som) . 
Likewise, for f € C 2m (SO(3)), the approximant converges in L^ with error 

Wf-TsfW^^Cp^WfWc^soo))- 
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Proof. For 1 < p < oo we can estimate the approximation error by 



\f-T s f\\ P 



< 



/ C m f(a)(k m (-,a) - 'Va(|,a)/e„ l (-,|))d/i(a) 
Jso(3) v |^ ' 

I |£ m /(a)| e fem (-,a)d//(a) 



< ||£? fcm ||p^ P ||/: m /|| P . 



We need to consider only two values of p, namely p = 1 and oo, because the 
others follow by interpolation. The L\ operator norm is controlled by M\ := 
su Paeso(3) / e * m a)d^i(s). Likewise, the norm is controlled by Mx, :— 
svlV x ^so{z) J ek m (x, cx)dfi(a). Since the estimate from Lemma [8] controls ek, m by a 
quantity involving uj(a~ 1 x) only, the error kernel is generated by a class function, 
and we can make a change of variable to see that the result is the same for L\ and 
ioo (and, hence, all other 1 < p < oo as well): 



max (Afi, Moo) < 



CP 2 



<SO{S) 

< c P 2m - 3 

< Cp 2 ' 



1 



o 

oo 



dist(-, a) 
P 

2m-A-L 



2m-4-L 



dp(a) 



sin — duj 



(1 + 9) 2m - 4 - L 2 d6<Cp 2 



Interpolating over 1 < p < oo we obtain 

H-Efem ||p->p < C p 2m , 
and, hence, ||/-T a /|| p < \\E km \\ p ^ p \\C m f\\ p < C p 2m \\C m f\\ 



□ 



Note 2. The size of the coefficients can be estimated as follows. Given a 
coefficient kernel a satisfying the CKCs with stability constant K, we have 



\A\\ ei{s) < K\\C m f\\ 



Indeed, 



E \ A *\ ^ f E |a(a,£)||£ m /(a)|d M (a) < W |£ m /(a)|d M (a). 

"'SO(3)^ e3 JSO(3) 

Applying Lemma [5] to the previous theorem gives the following corollary. 

Corollary 10. For centers H having density h < ho (with ho determined by L as 
in Lemma\^), there exists a coefficient kernel such that the approximant defined in 
(|24[) converges like 

\\f -T s f\\ p <Ch 2m \\f\\ w?Hsom] 
with coefficients satisfying the estimate 

\\A\\ ll{S) <2\\C m f\\ l . 
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5.2. Functions of Lower Smoothness. To treat more general functions, we make 
use of a common if -functional argument, one practically ubiquitous in approxima- 
tion theory. This will permit us to define approximants for functions in Besov 
spaces of lower order (and ultimately get good error bounds for such functions as 
well). 

Definition 11. For < s < 2m and 1 < p < oo, we define the Besov space 
^poo(^P)) as th e collection of functions in L P (S*0(3)) for which the following 
expression 

\\f\\ B ^:= sup t- s K(f,t) 
t>0 

is finite, where the K -functional K(f, ■) : (0, oo) — > (0, oo) is defined as 

K(f,t) :=inf {\\f-g\\ Lp+ t 2 ™\\g\\ wSm : S 6? p 2m }. 

For p = oo. the definition is the same after substituting L p (S*0(3)) by (7(5*0(3)) 
and Wl m by C 2m . 

There are abundant alternative characterizations of Besov spaces. We have cho- 
sen the preceding because it captures precisely the properties we need from a func- 
tion of lower smoothness. It is possible to follow the example of Definition [TJ by 
transporting the Besov norm from R 3 . For a general manifold this would not work, 
it is possible only because <S*0(3) is compact (specifically, it works because the 
diffeomorphisms and partitions of unity employed have finite cardinality- see |26[ 
1.11.3, Remark 2 and 2.4.7] for a discussion of this). We recommend 26, 1.11, and 
Chapter 7] and the references therein for background on Besov spaces on Riemann- 
ian manifolds. 

An immediate consequence of the definition is that, for an arbitrary function in 
and p > 0, there is g p € W p m simultaneously satisfying 

\\g p \\ w , m <2p s - 2m \\f\\ B s ^ (25) 
\\f-9 P \\L p <2p s \\f\\ Bioo . (26) 

In the setting of Theorem [9] (namely, with radius p), the approximant to / is then 
s/,s := T 3 g p . 

Note 3. It is important to point out that the function g p can be constructed by the 
same method employed in Section |2~31 to create the C°° approximant. Namely, one 
may employ a (7°° partition of unity (ipj)^^, a corresponding sequence of smooth 
diffeomorphisms, and a sequence of cut-offs (if: to generate an approximant of 
the form 

N 

where gj. p : K 3 — >• R approximates the function fj — (/ x ipj) o hj 1 £ B s p ^(IR 3 ) in 
the sense of (|2"5"|) and (f2"6"| . The Euclidean approximants g JiP can be constructed, 
for example, by the method [I] Theorem 5.33] or [6l Lemma 5.1]. 

Note 4. Again, we can estimate the size of the coefficients. For / g B s p x we have 

P||, l(3 )<2Xp s - 2m ||/|| B?i0o . 

This follows directly from Note [5] 
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Theorem 12. Let s < 2m and 1 < p < oo. Suppose the coefficient kernel a : 
S x SO(3) — » M satisfies CKCs with radius p, precision L > 2m, and stability 
constant K . Then for f € ^ there is Sf t s satisfying 

\\f-s f M\p<c P s \\f\\B^. 

Proof. This follows by applying Theorem |H] to g p and observing that (|23|) ensures 
II^IIp^P-^II/Hb;... " □ 

Applying Lemma [6] again, we obtain: 

Corollary 13. For centers S having density h < ho (with ho determined by L as 
in Lemma\B^, and for f € ^ f/iere is an approximant s/,s satisfying 

\\f-SfM\ P <Ch s \\f\\B koo 
with coefficients satisfying the estimate 

\\A\\ ei{s) <2h s - 2m \\f\\ B ^. 

6. Concluding Remarks 

The results of the last section have value beyond solving a challenging (if theo- 
retical) problem, and even beyond providing benchmarks for other approximation 
schemes. When approximation schemes based on projection (e.g., interpolation 
and least-squares projection, discussed in I4.2|) are stable, application of a Lebesgue 
lemma implies that such schemes are near best. Very recently, in |17j . it has been 
demonstrated that kernel interpolation on smooth, compact Riemannian manifolds 
is stable for a certain family of positive definite kernels. Under similar circum- 
stances in [16] , least-squares projection has been shown to be stable in L p for each 
1 < p < oo. In the setting of the 2-sphere, S 2 , it is known that these kernels 
provide precise (if theoretical) L p rates of approximation. Hence, interpolation on 
§ 2 (an easily implemented approximation scheme) and least-squares approximation 
inherit these precise rates for all levels of smoothness up to a saturation order of 
the kernel. 

The kernels employed in [17] [16] have variants on S*0(3) as well, and they are very 
similar to the kernels considered here. Both families are composed of fundamental 
solutions of elliptic differential operators which serve as reproducing kernels for the 
Sobolev spaces W™{SO(3)) (albeit for different, but equivalent inner products). 
This raises the strong likelihood that similar stability results hold for the kernels 
considered here, or that, by modifying the techniques of this article, precise error 
estimates can be shown for the kinds of kernels considered in jTTJ [TB] . The upshot 
is that precise error estimates for 'best' approximation may lead directly to precise 
error estimates for any stable linear approximation scheme obtained by projection. 

In addition to this, the approach we take has the potential to work in far greater 
generality: specifically on other structures like Lie groups and homogeneous spaces. 
It relies on two basic properties that, at first, seem to have little to do with groups, 
but rather with (compact) manifolds. These are: 

• the representation (1111) involving the fundamental solution for a differential 
operator: iterated, perturbed Laplace -Beltrami operators; 

• the ability to efficiently replace this fundamental solution by a linear com- 
bination of nearby translates incurring little error (Lemma [8|). 



SURFACE SPLINE APPROXIMATION ON SO(3) 



21 



Finding kernels that satisfy both of these simultaneously is a challenge; the second 
property seems especially difficult without some extra structure. However, when 
the manifold is a group or a homogeneous space, it is often possible to exploit 
some underlying symmetry. In this article we considered kernels, k(x,oc), derived 
from 'class functions', i.e. functions depending on a~ 1 x - meaning k(x,cx) = 
f(a~ 1 x) - but possessing an added symmetry; the function / depends only on the 
rotational angle of a~ 1 x (which on 50(3) is the same as the Riemannian distance). 
It is precisely this symmetry - ultimately embodied in the addition formula which 
relates the basic elements of harmonic analysis on the space (the eigenfunctions 
of the Laplace - Beltrami operator) to an orthonormal basis for the space of class 
functions, and which, in turn, can be related to certain orthogonal polynomials - 
that ultimately allowed us to effectively replace the kernel. 
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